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I . INTRODUCTION 


This document is the second annual report of the NASA 
Geodynamics Program. Its purpose is to inform interested 
agencies and individuals of the status, progress, and future 
plans of the program. The intention of the Geodynamics Program 
Office is to issue similar reports each year. 

The first annual report (NASA, 1980c) contained a section of 
background and historical information, to bring the reader up to 
1979 since the beginning of the program as the National Geodetic 
Satellite Program. The historical section is omitted here. For 
the first time we include a list of publications of research 
results supported by the Program (Appendix 3) . 


HIGHLIGHTS AND ACHIEVEMENTS OF 1980 

In 1980 the detailed planning and organization begun in 
1978-1979 was continued. Comments were received from recipients 
of the Geodynamics Program Plan (NASA, 1979a) and, informally, 
from several advisory committees: the National Academy of 

Sciences' Panel on Crustal Movements Measurements, Committees on 
Geodesy and on Seismology, Committee on Earth Sciences of the 
Space Science Board; and from NASA's Advisory Subcommittee on 
Geology and Geodynamics. 

Among the highlights of 1980 were the following: 

1. The Crustal Dynamics Project Plan was approved, and the 
Crustal Dynamics Project initiated (NASA, 1979c) . 

2. An agreement was signed between NASA, USGS, NSF, the 
Defense Mapping Agency, and the National Geodetic Survey (a NOAA 
agency) , establishing a coordinated Federal program for the 
application of space technology to geodynamics. The agreement 
provides for a Program Review Board at management level, and an 
Interagency Coordinating Committee at program office level. The 
ICC is a continuation of an informal coordinating committee which 
had been meeting since early 1979. 

3. The International Association of Geodesy, one of the 
constituent associations of the International Union of Geodesy 
and Geophysics, established a Commission on International 
Coordination of Space Techniques for Geodesy and Geodynamics. 

The Commission President is Prof. I. I. Mueller, a long-time 
participant in NASA's Geodynamics Program. Secretaries of CSTG 
are Prof. J. Kovalevsky (France), Prof. A. Massevitch (USSR), and 
Dr. E. A. Flinn, Chief Scientist of the NASA Geodynamics Program. 
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4. _ The Transportable Laser Ranging Station (TLRS) , designed 
and built by the University of Texas at Austin, was completed, 
checked out, and deployed to the Jet Propulsion Laboratory for 
TLRS/ARIES intercompar ison (Laser/VLBI Intercomparison Session 
IV) before being put into service in the Crustal Dynamics Proiect 
in early 1981. 

5. Comparison of VLBI and laser ranging baseline length 
determinations at observatory sites in New England, Texas, and 
California showed agreement between the two methods to better 
than approximately five centimeters. 

6. An Announcement of Opportunity for Crustal Dynamics Data 
Use was issued in October 1980. Approximately 90 proposals were 
received in response to this AO, including approximately sixteen 
foreign investigations. Investigations will be selected in 
mid-1981, and will begin in late 1981. 

7. Analysis of the laser ranging observations at Quincy and 
Otay Mountain (California) , a continuation of the previous San 
Andreas Fault Experiment measurements which have been acquired 
since 1972, showed good agreement with past measurements, and 
strengthen the conclusion that the baseline length between these 
two sites is decreasing at a rate of approximately 8 cm/year - a 
result consistent in direction but different in magnitude from 
that predicted by plate tectonics theory. 

The Magsat satellite re-entered the earth's atmosphere. 
The mission objectives were successfully accomplished. Magnetic 
field models were provided to USGS for use in preparation of 
magnetic charts and maps. 

9. The International Council of Scientific Unions (ICSU) 
approved a request from the International Union of Geodesy and 
Geophysics and the International Union of Geological Sciences to 
establish a successor project to the highly successful 
International Geodynamics Project, which had been approved for 
the period 1971-1980. The new project, called the International 
Lithosphere Program (ILP) , will be conducted through 1989 by the 
Inter-Union Commission on the Lithosphere. The scientific 
importance of space geodynamics is indicated by the fact that two 
of the five scientific objectives of ILP cannot be accomplished 
without a global program of space measurements of plate movement 
and deformation. Dr. R. A. Price of the Canadian Geological 
Survey is President of the Commission, and Dr. E. A. Flinn is 
Secretary-General . 

10. The Second Annual NASA Geodynamics Program Review was 
held in early 1980, and attended by most of the program 
participants. The five-day meeting consisted of presentations 
and discussion of scientific results. 
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11. A Gravsat User Working Group was formed by NASA, and met 
to continue the studies of a gravity field mission begun by the 
National Academy of Sciences in 1978. A Federal Gravsat panel, 
established under the auspices of the Satellite Geodesy 
Applications Board, developed a plan for joint agency 
participation in the Gravsat mission. 

12. A brassboard model of the Gravsat satellite- to-satellite 
tracking system and simulation of the drag-free system were 
initiated. (Recent results have confirmed that the required 
tracking accuracies are entirely feasible.) Computer simulations 
show that the Gravsat mission as presently configured can improve 
the accuracy of gravity field models by two to three orders of 
magnitude . 

13. NASA supported the 1980 Project MERIT observing campaign 
in August-October 1980 by concentrating VLBI and laser ranging 
observations in this period to the maximum extent possible. 

14. Fifty-nine independent baselines were measured in 1980 
using both laser ranging and VLBI; this number does not include 
baselines measured using European stations. For comparison, 55 
independent baselines were measured in 1979, and the Crustal 
Dynamics Project plans to measure 124 in 1981. 

15. An interagency plan for geodetic use of the DOD's Global 
Positioning System satellites was prepared (NASA et al., 1980a). 

16. A plan for laser development was published (NASA, 

1980a) . 
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II. GEODYNAMICS PROGRAM STRUCTURE AND STATUS 
A. PROGRAM SCOPE AND OBJECTIVES 


We summarize here the scope of the NASA Geodynamics Program 
and its objectives. For its historical development, see the 
Geodynamics Program Annual Report for 1979 (NASA, 1980c) . 

The objectives of the NASA Geodynamics Program are: (1) to 

contribute to understanding of the solid earth, in particular the 
long-term crustal processes associated with natural hazards and 
resources, and the structure and internal composition of the 
earth; (2) to facilitate the establishment of new geodynamics 
measurement services requiring precise position determination. 

Achievement of these objectives depends on the development, 
demonstration, and use of techniques for measuring crustal motion 
and deformation and on improved measurement of the earth's 
gravity and magnetic fields. The program includes research in 
modeling of crustal processes, earth rotational dynamics, and the 
geopotential field; this research is a necessary part of the 
program. Studies of new instrumentation and space mission 
concepts are also supported. 

The principal space techniques currently available for 
precise measurements of crustal motion and deformation, polar 
motion, and earth rotation are very long baseline interferometry 
(VLBI) using signals from extragalactic radio sources or from 
satellites, and laser ranging to the Moon and to man-made 
satellites such as Lageos. The VLBI and satellite laser ranging 
methods employ both fixed and mobile stations. 

The status of elements of the Geodynamics Program are 
described in the following sections. Further description of the 
research program results is given in Section III. 
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B. CRUSTAL DYNAMICS PROJECT 


Beginning in FY 1980, several smaller projects (SAFE, PPME, 
ARIES) were combined to form the Crustal Dynamics Project. This 
Project is managed by the Goddard Space Flight Center (GSFC) with 
the support of the Jet Propulsion Laboratory (JPL) . Other 
participants include the Smithsonian Astrophysical Observatory 
(SAO) , the University of Texas at Austin, the University of 
Hawaii, the Massachusetts Institute of Technology, and the 
Northeastern Radio Observatory Corporation (NEROC) . 

The objectives of the Crustal Dynamics Project are to measure 
and model: 

1. Regional deformation and strain changes related to 
earthquakes at the plate boundary in the western U.S. and Alaska; 

2. The present relative velocity of major plates, with 
emphasis on the North American and Pacific Plates and those 
plates interacting with the Pacific Plate; 

3. Internal deformation of the North American and Pacific 
Plates; 

4. Regional deformation in regions whose tectonic setting 
is similar to that of the Western U.S. and Alaska. 


The Crustal Dynamics Project is making use of existing 
systems (developed to support SAFE, PPME and ARIES) , upgrading 
the performance of those systems, and developing new systems in 
order to conduct a systematic program of crustal movement and 
deformation measurements through 1986. This NASA Project is an 
integral part of the multi-agency Federal program for the 
application of space technology to crustal dynamics and 
earthquake research. Cooperative arrangements with other 
countries either having or planning similar capabilities has 
resulted in extension of project activities to a broad-based 
global research program in geodynamics. Investigators selected 
under the October 1980 Announcement of Opportunity for Crustal 
Dynamics will form an Investigator's Working Group and will 
become part of the Project beginning in 1982. 

The location and dates of operation of funded laser and VLBI 
systems and cooperating stations in other countries is given in 
Table II-l. 
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TABLE II-l 


STATIONS AND LOCATIONS FROM WHICH LASER AND/OR VLBI DATA 

ARE AVAILABLE * 


Year Geodetic East 


74 

75 

76 

77 

78 

79 

80 

station 

Latitude 

Longitude 








Goddard Laser Data 





X 

X 

X 

X 

X 

X 

X 

Greenbelt (GSFC) Md, 

39‘ 

^01' 

- 76‘ 

^50' 

X 


X 



X 


Quincy, California 

39 

58 

-120 

56 

X 


X 

X 

X 

X 


Otay Mt . , Cal . 

32 

36 

-116 

50 


X 

X 


X 



Bermuda 

32 

21 

- 64 

39 


X 

X 


X 



Grand Turk 

21 

28 

- 71 

08 


X 

X 

X 

X 

X 

X 

Patrick AFB, Fla. 

28 

14 

- 80 

36 



X 



X 


Bear Lake, Utah 

41 

56 

-111 

25 




X 

X 

X 

X 

Westford, Mass 

42 

37 

- 71 

29 








(Haystack Obs.) 








X 

X 

X 

X 

Big Pine, Calif. 

37 

14 

-118 

18 








(Owens Valley RO) 









X 

X 


Goldstone (DSS-14) , Ca 

. 35 

25 

-116 

53 






X 

X 

Goldstone (DSS-13) , Ca 

. 35 

15 

-116 

48 






X 

X 

American Samoa 

-14 

20 

-170 

44 






X 

X 

Yaragadee, Western 












Australia 

-29 

03 

115 

21 






X 

X 

Kwajalein 

9 

24 

167 

29 






X 

X 

Ft. Davis, Texas 












(McDonald Obs.) 

30 

41 

-104 

01 







X 

Haleakala, Maui, 












Hawaii (LURE Obs.) 

20 

42 

-156 

15 








SAO Laser Data 






X 

X 

X 

X 

X 

X 

Mt. Hopkins, Arizona 

31 

41 

-110 

53 


X 

X 

X 

X 

X 

X 

Arequipa, Peru 

-16 

28 

- 71 

30 


X 

X 

X 

X 

X 

X 

Natal, Brazil 

- 5 

56 

- 35 

10 



X 

X 

X 

X 

X 

Orroral Valley, 












Australia 

-35 

38 

148 

57 








European Laser Data 








X 

X 

X 

X 

Kootwijk, Netherlands 

52 

11 

5 

49 




X 

X 

X 

X 

Wettzell, W. Germany 

49 

09 

12 

53 




X 

X 

X 


San Fernando, Spain 

36 

28 

- 6 

12 








Lunar Laser Data 





X 

X 

X 

X 

X 

X 

X 

Ft. Davis, Texas 












(McDonald Obs.) 

30 

41 

-104 

01 



X 





Haleakala, Maui, 












Hawaii (LURE Obs.) 

20 

42 

-156 

15 






X 

X 

Orroral Valley, 












Australia 

-35 

38 

148 

57 
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TABLE II-l (continued) 


Year Geodetic East 


74 

75 

76 

77 

78 

79 

80 

station Latitude 

Longitude 








Observatory VLBI Data 





X 

X 

X 

X 

X 

X 

X 

Westford, Mass. 












(Haystack Obs.) 

42 

37 

- 71 

29 

X 

X 

X 





Goldstone (DSS-14 ) , Ca. 

35 

26 

-116 

53 

X 

X 

X 

X 

X 

X 

X 

Greenbank, W. Va. 












(NRAO) 

38 

26 

- 79 

50 

X 



X 

X 

X 

X 

Onsala Space Obs., 












Onsala, Sweden 

57 

24 

11 

56 



X 

X 

X 

X 

X 

Big Pine, Calif. 












(Owens Valley RO) 

37 

14 

-118 

17 






X 

X 

Efflesberg, W. Germany 

50 

31 

6 

53 







X 

Ft. Davis, Texas 












(Harvard RAS) 

30 

38 

-103 

57 








Mobile Station VLBI Data 





X 

X 

X 

X 


X 


Goldstone (DSS-14) , Ca. 

35 

26 

-116 

53 




X 


X 

X 

Goldstone (DSS-13) , Ca. 

35 

15 

-116 

48 

X 

X 

X 

X 

X 

X 

X 

Pasadena, California 












(Jet Prop. Lab. ) 

34 

12 

-118 

10 


X 

X 

X 

X 

X 

X 

Big Pine, California 












(Owens Valley RO) 

37 

14 

-118 

17 


X 


X 




Malibu, California 

34 

04 

-118 

39 



X 

X 



X 

Palos Verdes, Calif. 

33 

45 

-118 

24 



X 

X 




Pearblossom, Calif. 

34 

31 

-117 

55 




X 



X 

La Jolla, California 

32 

52 

-117 

15 




X 




San Francisco, Calif. 

37 

48 

-122 

27 






X 


Quincy, California 

39 

58 

-120 

56 






X 


Otay Mt . , Calif . 

32 

36 

-116 

50 


* Some data taken earlier than 1974 is also available. 
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In 1980 an engineering version of a 4m mobile VLBI was 
refurbished for field use and procurement was initiated (for 
delivery in early 1983) of the first mobile VLBI unit 
specifically designed for field operations (Figure 1) . 

Procurement was also initiated for a transportable VLBI data 
system which can be used to directly convert radio astronomy 
antennas to geodetic VLBI facilities. Also in 1980, the first of 
an improved version of mobile lasers (Transportable Laser Ranging 
Station - TLRS-1; see Figure 2) was completed and tested and 
integration of a modularized laser ranging system (TLRS-2; see 
Figure 3) was initiated. 

The ranging precision of the NASA lasers at present is 
typically better than 10 cm: the NASA Stalas system has a 
precision of 1-2 cm, three of the Moblas have precisions of 3-5 
cm, five of the Moblas systems have precisions of 6-9 cm, and the 
TLRS have precisions of 1-2 cm. In 1980 an interim modification 
the eight Moblas was initiated to improve ranging 
precision to the 2-5 cm level. This modification will be 
completed in 1981. 

Discussions were also initiated with the Australian Division 
of National Mapping (Natmap) for modification of the Natmap Lunar 
Laser Ranging Facility at Orroral for satellite ranging. 

The development of laser ranging systems for crustal dynamics 
applications is desctibed in a plan issued in December 1980 
(NASA, 1980b) . 

The current precision of VLBI systems for position 
determination is 3 to 6 cm. Upgrading of the VLBI technique is 
planned to achieve precisions of 1 to 3 cm. A key factor in 
achieving this performance is accurate correction for path length 
changes caused by water vapor in the atmosphere. Water vapor 
radiometers capable of reducing this error source to about 2 cm 
have been developed. Studies of improvements in instrument 
calibration and atmospheric modeling designed to reduce the error 
to the 1 cm level were initiated in 1980. 

NASA plans for development of the VLBI method for geodetic 
applications and interplanetary spacecraft tracking are outlined 
in a document issued in December 1979 (NASA, 1979b) . 

In mid-1980, a five-station experiment (Figure 4) was 
conducted with Moblas systems at Haystack, Massachusetts; Ft. 
Davis, Texas; Goldstone, California; Owens Valley, California; 
Stalas at GSFC; and VLBI facilities at Haystack, NRAO, West 
Virginia; Ft. Davis, Goldstone, and Owens Valley. Doppler data 
were acquired by NOAA at each of these sites and at five other 
locations in the U.S. The goal was to achieve agreement between 
the VLBI and laser systems at the 5 cm level or better. 
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In October 1980 the OVRO-Goldstone-JPL baselines in 
California were measured using TLRS-1, the 4 meter mobile VLSI at 
JPL, and Moblas stations at the OVRO and Goldstone VLBI sites. 

The objective is to verify the accuracy of the mobile VLBI 
system. 

In 1979 and for part of 1980, laser ranging provided initial 
baselines, with which later measurements will be used to 
determine the relative motion of the North American and Pacific 
plates. The sites were located in Hawaii, American Samoa, and 
Kwajalein in the Pacific; and Ft. Davis, Texas, Goddard Space 
Flight Center, Maryland, and Owens Valley, California, in North 
America. 

Measurements of the relative plate motion across the North 
Atlantic using VLBI techniques also began in 1979. Four sites in 
North America (Westford, Massachusetts; Greenbank, West Virginia; 
Owens Valley, California; and Ft. Davis, Texas) have been tied to 
European stations (Onsala, Sweden and Bonn, West Germany) . 
Currently, measurement of intersite distances are being obtained 
at about one-month intervals between the Ft. Davis and Haystack 
stations on the North American Plate and the Onsala station on 
the Eurasian Plate, as part of the NOAA Polaris tests. VLBI 
measurements have already been carried out between Haystack, 
Greenbank, Owens Valley, and Goldstone, and these place 
constraints of less than one cm/year on the internal deformation 
of this part of the North American Plate (Figure 12) . 

For regional deformation and North American Plate deformation 
studies, about 80 sites throughout the continental United States 
have been recommended by representatives of USGS, NGS, NASA and 
interested university and government groups. More than 59 of 
these locations are in the western United States (Figure 5) and 
35 of these in California, as shown in Figure 6. These 
recommendations have been reviewed in terms of scientific value, 
priorities, implementation feasibility, availability of 
measurement systems, and measurement strategy, in order to 
establish the specific sites to be used. Not all the sites will 
be occupied with the same frequency, since the tectonic problems 
themselves are varied and the expected rates of motion differ 
from one region to another. 

Both laser ranging and VLBI measurements have been made over 
the past nine years at several locations in California. These 
sites have been incorporated into the regional deformation 
measurement network, in order to make use of the past 
measurements. Included are eight mobile VLBI sites, five mobile 
laser sites, and two VLBI base stations. In addition to these 
locations, other sites in or near California will be visited on a 
regular basis by mobile systems. These sites were selected 
partially to link together existing geodetic networks operated by 
NGS and USGS. 
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The U.S. regional deformation measurements, which will begin 
in 1981, will involve about 25 site visits per year. The U.S. 
sites will evolve into a National Crustal Motion Network which 
will be monitored by NOAA, beginning in FY 1984. 
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C. LASER NETWORK OPERATIONS 

The NASA satellite laser network^ which supports the Crustal 
Dynamics Project, consists of the eight Moblas operated by GSFC 
and four fixed lasers operated by the SAG. A high accuracy fixed 
laser (Stalas) is located at GSFC. In 1980, arrangements with 
the USAF for operation of a laser in Florida (Ramlas) were 
discontinued. During the year Moblas lasers were located at 
GSFC, Haystack, Massachusetts; Ft. Davis, Texas; Owens Valley, 
California; Goldstone, California; Kwajalein Island, American 
Samoa; and Yarragadee, Australia. A Moblas was also collocated 
with the fixed laser at the LURE Observatory at Haleakala, Maui, 
Hawaii, and remained there while the Hawaii laser operations were 
interrupted for refurbishment and upgrade of the facility. 

SAO lasers were operated in Arequipa, Peru; Natal, Brazil; 
and Canberra, Australia. A fourth SAO laser at Mt. Hopkins, 
Arizona, was used to test planned improvements to SAO lasers. 
These improvements which will be incorporated by earlier 1982 are 
expected to improve ranging performance to Lageos from 10 cm to 
3-5 cm. 


D. RESEARCH AND TECHNIQUE DEVELOPMENT (RTD) 


The RTD portion of the NASA Geodynamics Program for 1980 
included five subprogram elements. These are: Global Earth 

Structure and Dynamics; Lithospheric Structure and Evolution; 
Regional Crustal Deformation Modeling; Geopotential Field 
Modeling; and Advance Systems and Mission Studies. Portions of 
these activities are carried out through internal research at 
NASA Centers (GSFC, JPL, and MSFC) and through organizations 
selected on the basis of proposals submitted in response to an 
Applications Notice for Geodynamics. For FY 1980, 25 of the 78 
proposals received in response to AN;OSTA:80-B were selected for 
funding. Results of research funded under Geodynamics ATD are 
described in Section III. 

The Advanced System and Mission Studies included continuation 
of development of SERIES, a geodetic receiver concept using 
telemetry from the Global Positioning Satellite; a cryogenic 
gravity gradiometer, and a shuttle time and frequency transfer 
experiment. New activities were initiated for studies of 
atmospheric effects on laser and microv;ave systems, and an 
airborne laser ranging system for geodetic measurements. 
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E. MAGNETIC FIELD SATELLITE (MAGSAT) 

Magsat (Figure 7) was launched on October 30/ 1979, and 
re-entered the atmosphere on June 11, 1980. It carried a three 
axis vector magnetometer and a scalar magnetometer. The 
objectives of the mission were: 

1. To obtain an accurate description of the earth's time 
varying main magnetic field for use in updating and refining 
world and regional magnetic charts and maps. 

2. To compile global crustal magnetic anomaly maps and to 
interpret them in terms of geological/geophysical models of the 
crust for natural resource assessment purposes. 

During its seven and one-half months in orbit, Magsat 
provided the most accurate measurements of the global field ever 
obtained and the first measurements of the vector field in low 
earth orbit. Initial field models were provided in 1980 to 
Magsat investigators. Final processing of all Magsat data is 
expected to be completed in 1982. 


F. EARTH GRAVITY FIELD SURVEY MISSION ( GRAYS AT) 

Requirements for a gravity field mission were studied by the 
Committee on Geodesy of the National Academy of Sciences in 1979 
(NAS, 1979) . In 1980 a Gravsat User Working Group (GUWG) was 
formed by NASA to continue the studies conducted by the NAS and 
to address questions relating to the feasibility of the mission 
and the recovery of gravity information (NASA, 1980b) . In 
addition, a Federal Gravsat Panel was initiated under the 
auspices of the Satellite Geodesy Applications Board to develop a 
plan (NASA et al., 1980b) for joint agency participation in the 
Gravsat mission. 

The Gravsat mission will use satellite-to-satellite tracking 
between two low-altitude (160 km) drag-free satellites in polar 
orbit (Figure 8) . A two-way coherent microwave Doppler system 
will be used to achieve tracking accuracies of one micrometer per 
second. A brassboard of this tracking system and 
six-degree-of-f reedom simulations of the drag-free system were 
initiated in late 1980. These studies are planned to be 
completed in mid 1981. 
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G. INTERAGENCY COORDINATION 

The NASA Geodynamics Program activity was designed to fit the 
interests and objectives of other Federal agencies, and its 
activities are closely coordinated with those of other agencies. 
In particular, it is an objective of the NASA program to assist 
other Federal agencies in the establishment of improved or new 
geodetic services using space technology. 

In September 1980 an interagency agreement was consummated by 
NASA, NOAA, USGS , DMA, and NSF , providing for development of 
space technology for crustal dynamics and earthquake research. 

The agreement establishes a high-level Program Review Board and 
an Interagency Coordinating Committee, and requires development 
of a Federal Implementation Plan for the joint development and 
use of space systems for the acquisition and analysis of data, 
for the transition of new technologies from R&D to operational 
status, and for expansion of the national effort to a global 
program of geodynamics research. Specifically, three new 
services are planned; (1) a National Crustal Motion Network by 
NOAA beginning in 1984; (2) an improved Polar Motion and Earth 
Rotation Service (Polaris) jointly developed by NASA and NOAA and 
operational in 1983; and (3) a Local Crustal Motion Network, 
based initially on the use of SERIES, operated by the USGS 
beginning in 1983. 

In view of the potential importance of GPS to civilian 
satellite geodesy, an interagency group including NOAA, DOD, 

NASA, and USGS has prepared a plan for the development and 
testing of GPS receiver concepts (NASA et al., 1979). This plan 
will coordinate the activities of the several Federal agencies 
involved in developing GPS applications, and will provide a basis 
for^ testing these concepts and eventually selecting the optimum 
method, based on costs and performance, for general use. As 
outlined in the plan, NASA, in proceeding with the development of 
SERIES (which does not require knowledge of the signal code) , 

NOAA, DMA, and USGS are supporting development of a "tri-agency" 
receiver based on reconstruction of the GPS carrier phase. 
Engineering prototypes of SERIES have been built and will be 
tested by mid-1981. Prototypes of the "tri-agency" receiver are 
expected to be available in late 1983. A decision as to which 
receiver will be procured by the participating agencies will be 
made in 1984, and should result in fieldable units in 1985. 
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H. INTERNATIONAL COORDINATION 

1. Space Systems . In 1980, several laser systems operated 
by other countries were cooperating with the NASA program. These 
include a high-accuracy lunar/satellite laser in the Federal 
Republic of Germany and satellite lasers in France, The 
Netherlands, Spain, Greece, and Egypt. England, Austria, Japan, 
and Italy are either planning satellite lasers or have them under 
development. The Lunar Laser Ranging Station at Orroral Valley, 
Australia, is also a cooperating station. 

Cooperating VLBI systems operated by other countries are 
those in Onsala, Sweden, and Wettzell, Federal Republic of 
Germany. A VLBI system under development in Japan is expected to 
be operational in 1983, to participate in joint US-Japan 
experiments in 1983-84. The Federal Republic of Germany is 
planning to procure a VLBI facility dedicated to geodetic use, 
and other countries (e.g. France, Italy, and England) are 
considering development of observatory VLBI facilities. 

2. Project Merit . Commissions 19 and 31 of the 

International Astronomical Union have developed plans for a 
special period of international collaboration in the Monitoring 
of Earth Rotation and Intercomparison of the Techniques of 
observation and analysis. Project Merit. The objectives of the 
projects are: (1) to foster the development of new techniques 

for the measurement of the variations in the rate and axis of 
rotation of the earth, (2) to obtain precise data on earth 
rotation in order to increase our understanding of the causes and 
effects of the variations, and (3) to make recommendations on the 
observational basis and organizational arrangements for future 
international services on earth rotation. 

An initial period of observations to test techniques and 
arrangements for international cooperation occurred from 
September 1 to November 1, 1980. During this period, special 
efforts were made by NASA and other organizations to acquire 
laser and VLBI observations of polar motion and earth rotation 
and to accelerate the processing of data. Data analysis is now 
under way, and results will be reported in due course. 
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I. GEODYNAMICS PROGRAM FUNDING 

Total NASA funding for support of the Geodynaniics Program is 
$25. 3M in FY 1980 and $26. 9M in FY 1981 (Table II-2) . In 
general, about 40% of this funding is in support of regional 
crustal deformation studies in the U.S. Studies of plate motion 
and deformation are approximately equally supported and account 
for about 50% of the total. The remainder (10%) supports 
measurement and modeling of the geopotential fields. A 
significant portion of Program funds (approximately 70%) is 
devoted to system development and operations, largely through 
industry groups and NASA Centers. Research supported at 
universities, NASA Centers, and other organizations represents 
about 30% of the total. 


TABLE I I- 2 


NASA GEODYNAMICS PROGRAM FUNDING SUPPORT 
(IN MILLIONS) 



FY 1980 

FY 1981 

Crustal Dynamics Project 

$9.6M 

$11. 6M 

Laser Network Operations 

10.8 

10.9 

Research and Technique Development 

3.0 

3.2 

Magsat Mission 

1.6 

0.5 

Gravsat Studies 

0.3 

0.7 


$ 25. 3M 

$ 26. 9M 
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III. RESEARCH STATUS AND RESULTS 


Five main areas of research are conducted under the 
Geodynamics Research and Technique Development (RTD) Program. 
Most of the research results described this section were 
presented at the Third Annual Geodynamics Program Review, held 
the week of January 26, 1981, at Goddard Space Flight Center. 


A. PROGRAM ELEMENTS AND OBJECTIVES 


Before discussing the individual research results, we 
describe the objective of each of the five main elements that 
make up the Geodynamics RTD Program. The elements of Geodynamics 
research are: 

1. Global Earth Structure and Dynamics . The objective is to 
improve our understanding of the dynamics of the earth by 
development of models of polar motion and earth rotation, global 
plate motion, and the dynamics of the core; and to improve our 
understanding of the global structure of the earth, including its 
crustal magnetization, gravity field, and the evolution of the 
crust and lithosphere. Formulation of a standard dynamic earth 
model will be attempted under this program element. 

2. Regional Crustal Deformation Modeling . The objective of 
this element is to conduct modeling studies of crustal 
deformation in various tectonic settings. These studies are 
needed to assist in determining what measurements are required 
and to provide a proper perspective for analysis of data acquired 
by the Crustal Dynamics Project. The main types of models 
involved are earthquake mechanism, accumulation of stress and 
strain, and vertical motion. 

3. Lithosphere Structure and Evolution . In this element 
studies are conducted which are related to the processes which 
have formed the lithosphere, in order to gain a better 
understanding of the dynamic mechanisms that are active at 
present. Under this program element models are being developed 
that relate to the current state, origin, evolution, and dynamics 
of the lithosphere. The objective is to derive plate-driving 
mechanisms, including mantle convection, to explain current plate 
motion measurements. Emphasis is placed on understanding the 
nature of subduction and collision zones and lithosphere/mantle 
rheology. 

4. Geopotential Field Models . The objectives of this 
element are to develop gravity and magnetic field models, to 
investigate data analysis techniques and software systems, and to 
support the Non-Renewable Resources and Ocean Research Programs 
as well as other elements of the Geodynamics Program. Emphasis 
is placed on use of satellite altimetry and determination of the 
maximum resolution that can be achieved in the gravity field 
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analysis with existing data. Magnetic field models (including 
secular variation) are developed in this area of the program. 

The development of a data base for ancillary information required 
for geodynamics investigations, such as gravity anomalies, 
topography, and bathymetry, is maintained under this element. 

5. Advanced Geodynamics Instruments . These studies support 
the development of new methods for making geodynamics 
measurements using space systems and techniques. 


B . RESEARCH RESULTS 


Although there is some overlap between these various elements 
of the Geodynamics RTD Program, they form a logical framework for 
this research. Appendix 1 gives a listing of the current 
investigators, their affiliations, and the title of their current 
investigations. 

1. Global Earth Structure and Dynamics . 

This topic falls naturally into two categories: polar motion 

and earth rotation, and internal dynamics of the earth. We 
describe current results in rotational dynamics first, and then 
internal earth structure. 


a. Polar Motion and Earth Rotation 

One of the significant results comes from the analysis of 
over four years of upper-air zonal wind data to compute the 
angular momentum of the global atmosphere, and the resulting 
changes in length of day. In addition to distinct seasonal and 
interannual differences, higher frequency fluctuations are 
evident. Independent measurements of length of day determined by 
BIH, Lageos, and lunar laser ranging appear to correlate well 
with the atmospheric data (Figure 9) . The atmospheric data will 
be further studied to determine the regions that contribute most 
to the observed changes in the total angular momentum. 

Several groups are studying laser ranging data to determine 
polar motion and variations in the earth's rotation rate. Both 
long-arc and short-arc solutions for Lageos have been determined, 
and the results compared with BIH and Doppler solutions. The 
satellite laser ranging data compare favorably with the other 
solutions, particularly during the MERIT short campaign in late 
1980: differences of 0.02 arcseconds and 0.3 milliseconds in 

length of day are typical. 
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McDonald Observatory is still the only routinely operating 
lunar laser ranging observatory, so the types of analysis that 
can be performed with this data are limited. Variation of 
latitude and UTO have been determined from the McDonald lunar 
laser ranging observations for the period October 1970 to October 
1979, and these estimates have been compared with those 
determined from the pole position and UTl values obtained using 
other techniques. The McDonald results agree to within 5.2 
milli-arcseconds with the Doppler pole positions and to within 
1.6 milliseconds with the BIH UTO. The lunar laser ranging UTO 
values were used in conjunction with BIH values of polar motion 
to determine values of length of day and have been compared with 
changes in length of day inferred from a time series of global 
atmospheric angular momentum; the approximate 50 day oscillation 
recently discovered by Feissel and Gambis at the BIH have appear 
in both data types and agree closely in both amplitude and phase. 

Results of analysis of McDonald Observatory lunar laser 
ranging data include a value of GM of 398600.45+0.02 km^/sec^ and 
a secular acceleration of the lunar orbital mean longitude of n = 
-23.6+1.5"/century^, which yields a Q of 12.4 at semidiurnal 
frequencies. Tidally driven terms in UTl have coefficients k/c 
of 0.98+0.18 and 0.92+0.12 at fortnightly and monthly periods, 
respectively, compared with the theoretical value of 0.94+0.04. 
Also resulting from the solution are geocentric coordinates of 
McDonald accurate to 30 cm, information on the lunar orbit, lunar 
rotation, and reflector coordinates. 

Earth orientation parameters derived from Lageos observations 
were investigated in conjunction with data available from the 
Green Bank connected element interferometer, Doppler satellite 
tracking, and the classical optical instruments of the U.S. Naval 
Observatory. The results indicate that data contributed from 
each of the above techniques may be useful in a routine service 
for users. All sets of data may be included in the U.S. Naval 
Observatory algorithm providing earth orientation parameters and 
their predictions on a weekly basis. 

The comparison of determinations of UTl-UTC by independent 
techniques in 1979 led to the detection of fluctuations in the 
length of day with a total amplitude of 0.35 milliseconds and a 
period of 55 days. As a consequence, the estimated uncertainties 
(about 1 millisecond) of the five-day values of UTl derived by 
the BIH from classical astrometry must be considered realistic. 
The analysis of this series since 1967 shows that transient 
fluctuations are often present in the earth's rotation; their 
amplitudes are several milliseconds and periods 50 to 80 days. 

Using VLBI data from 1972 through July 1980, UTl and the 
x-component of polar motion have been determined with formal 
(one-sigma) errors since 1976 of typically 0.2 milliseconds for 
UTl and 3 milli-arcseconds for the x-component. Preliminary 
individual one-day solutions from North American data for the 
MERIT campaign (Figure 10a and 10b) show smooth changes of earth 
orientation with peak-to-peak variations of 1.3 milliseconds in 
UTl and 15 milli-arcseconds in the x-component (Figure 10a) over 
an interval of one week. 
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An accurate description of the elastic-gravitational 
deformation of a rotating, slightly elliptical earth subject to 
external gravitational forces has been developed. Some immediate 
consequences are a more precise theoretical nutation series and a 
uniformly valid description of the spatially varying Love numbers 
needed to characterize the behavior of an earth which is not 
spherically symmetric. 

Theoretical estimates of the Chandler period, accounting for 
the earth's slightly non-hydrostatic equilibrium shape, have been 
refined and have produced the first correct estimate of Chandler 
damping caused by anelasticity in the mantle and crust. These 
results served to significantly tighten bounds on the allowable 
frequency dependence of anelasticity over the period range of one 
hour to fourteen months. 

The complex interaction of pressure, gravity, and the 
magnetic field in the core may influence the dynamical motion of 
the earth. Core effects have been taken into account in 
computing corrections in the amplitude of the 18.6 year nutation 
obliquity of the earth. 


b. Internal Dynamics 

Results of numerical investigations of convection in a plane 
layer with temperature and depth dependent viscosity over a large 
parameter range show that when both top and bottom boundaries are 
stress free, the cold high-viscosity top thermal boundary layer 
has a dramatic effect on the cell aspect ratio and vertical 
convective flow structure. The precise manner in which the 
shallow high-viscosity layer is stripped off may critically 
affect model results, so that viscosity variations across the top 
thermal boundary layer may result in significant differences 
between the interior temperature, surface deformation, and 
gravity anomalies associated with constant and variable viscosity 
convection. 

Several other models being studied include variable rheology 
with depth. Using these models, it was found that a shallow 
region of low viscosity beneath a fault can substantially modify 
the predicted post-seismic vertical displacements and can lead to 
a large buildup with time in the post-seismic shear stress ahead 
of the fault tip. The uplift data following the 1964 Alaskan 
earthquake are consistent with models including both the 
subducted slab and a wedge of low viscosity material above the 
slab. Three-dimensional viscoelastic models of deformation 
following large thrust earthquakes indicate that surface 
horizontal velocities of 1 cm/yr or greater may be measurable 
within 1000 km of the plate boundary. 
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Locations have been determined of all the laser observatories 
which were operating during the period 1976 through June 1979. 

The baselines between the Greenbelt, San Diego and Areguipa 
(Peru) stations show the best repeatability and these stations 
are believed to have been the most consistent and accurate 
systems during these periods. Baseline consistencies between 
these stations are at the 5 cm level. The heights estimated in 
these solutions show consistencies at the 10 to 15 cm level. 


2 . Regional Crustal Deformation Modeling 

Several themes characterize current work on crustal movements 
and crustal deformation. One is the growing need for extension 
and refinement of the body of concepts that comprise the plate 
tectonics hypothesis. Near active plate margins the complexities 
of the real earth must be taken into account in constructing 
models, since regional-scale deformation of the plates is known 
to be occurring — for example, in California, where the global 
model of Minster and Jordan predicts 5-6 cm/year movement between 
the Pacific and North American plates, in contrast to 
approximately 3 cm/year of creep in Central California, no 
observed movement in Southern California, and approximately 8 
cm/year observed point-to-point between Quincy in the Northern 
Sierra and Otay Mountain on the Pacific plate near San Diego 
(Figure 11) . 

Another theme is the need for detailed study and modeling of 
the rheology of the lithosphere and asthenosphere in order to 
explain regional deformation. Finally, it is clear that in order 
to make progress in understanding the important geodynamical 
problems, the multiple disciplines of geodesy, geophysics, and 
geology will all have to be involved. Theoretical models draw on 
concepts and data from field investigations, laboratory analysis, 
and numerical experimentation. 

Kinematic models of plate motion continue to be refined; 
sources of complexity now being considered include possible 
non-rigid behavior (e.g., the Indian-Australian plate, which is 
being sheared along the Ninetyeast Ridge) , inclusion of 
small-area oceanic plates, and departures from ideal behavior 
near active margins (e .g . , inclusion of microplates at the north 
coast of South America and the west coast of North America in 
order to explain observed data) . 

Repeated observations of baseline lengths between VLBI 
observatories in the US and Europe show no changes at the level 
of 1 cm/year over the 4000 km baseline between Massachusetts and 
California since 1972 (Figure 12) and no change at the same level 
between West Virginia and Massachusetts. Predicted movement 
between North America and Europe is sufficiently slow that 
observations up to the present are simply initial-epoch 
measurements . 
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Vertical tectonic movements in presumably stable areas 
continue to be the subject of study. Comparisons between 
geological and geodetic evidence for movement indicate episodic 
behavior in the long term. Previous tectonic interpretations of 
releveling differences may be in error, however, since they have 
not always taken into account non-tectonic effects such as water 
withdrawal and systematic errors. 

Regional geologists continue to demonstrate the fruitfulness 
of collaboration between their discipline and those of geodesy 
and geophysics; examples are studies now being undertaken in the 
Alps, the Caribbean, China, and Asia Minor. 

Several studies continue of time-dependent regional surface 
deformation associated with earthquakes and fault motion. Creep 
in both the lower lithosphere and asthenosphere is a possible 
mechanism for significant postseismic surface motion and 
restressing of the earthquake slip zone. It appears that 
significant postseismic motion is possible at distances up to 
several hundred kilometers from a major fault, on time scales up 
to a few centuries - depending on the viscosity profile within 
the earth. One of the most important aspects of regional crustal 
deformation measurements is the light they may shed on the 
rheological properties of the asthenosphere. 

A new model for coseismic strain release includes the 
existence of a thin viscous decoupling layer at intermediate 
depths in the crust; preliminary results presented at the program 
review appear interesting. 

Seismological monitoring in the vicinity of McDonald 
Observatory in southwestern Texas has revealed active tectonic 
features which are not indicated by topological trends in the 
area. Studies of tectonics in several areas, using orbital 
photographs and Landsat images, are being conducted by several 
groups. 


3 . Lithospheric Structure and Evolution 

Several studies are being conducted of the relationship 
between regional geoids and earth structure, particularly in 
subduction zones. In the Aleutian/Alaskan region, models 
indicate that gravitationally apparent mass excesses do exist in 
subducting slabs, and are less than those generally predicted by 
thermal conduction models. 

Residual geoid .anomaly profiles for island arcs are being 
used to constrain numerical models of flow and temperature in 
subduction zones with a view to obtaining a better understanding 
of the deep structure of these regions and to address the 
question of the depth to which the return flow extends. Initial 
results obtained using a simple Newtonian model show a broad zone 
of downwelling centered on the trench and predict geoid anomalies 
of an order of magnitude larger than those actually observed. 



22 


Residual geoid anomaly maps show a positive anomaly 
associated with residual depths south of Australia, and a 
negative anomaly associated with the West Philippine Basin. 

These residual geoid anomalies are consistent with anomalous 
surface wave velocities and with a thermal origin, but not with 
residual depth. If both residual depth anomalies are caused by 
the same mechanism, then compositional variations in the 
uppermost mantle, rather than thermal variations, are probably 
the cause. 

Geoid lows and anomalously shallow sea floor occurring behind 
drifting continents (south of India, southwest of Australia, and 
east of North and South America) may be continental wake 
phenomena. Several differences occur between the oceanic and 
continental upper mantle (viscosity, composition and lithospheric 
thickness) that may be important for creating continental wakes; 
continental tilting and shedding of high-density mantle are 
probably not important. Less dense mantle rising in the wake of 
the continent would produce geoid and topographic highs if the 
mantle beneath the oceanic lithosphere behind the continent were 
of uniform viscosity. However, if a zone of low viscosity 
beneath the oceanic plate reduces the normal stress on the base 
of the lithosphere and the resulting uncompensated topography, a 
geoid low may occur in association with a topographic high, as 
observed in continental wakes. 

The global intraplate stress field predicted by plate-driving 
force models has been calculated as a means of testing such 
models against available observations on midplate stress 
orientations. The typical magnitude of intraplate stresses is 
several hundred bars, and local sources of stress of larger 
magnitude can dominate regional stress. The average changes in 
stress along faults during earthquakes are typically a small 
fraction of the total stress. For regional-scale problems, a 
linear viscoelastic but spatially variable rheology for the 
lithosphere and asthenosphere has been used to explore the 
effects of lateral variation in viscosity for both strike-slip 
and dip-slip earthquakes. 

Several investigators are combining data types for 
multi-disciplinary studies of particular regions; for example, 
seismic, magnetic, and gravity are being combined with geological 
and geodetic data for studies in North America. 


4. Geopotential Field Models. 

Significant progress has been made in modification and 
improvement of the Geodyn Orbit Determination and Parameter 
Estimation Program System which is used to recover geodetic and 
geophysical parameters from satellite data in a state-of-the-art 
manner. For example, time-dependent drag coefficients and 
altimeter crossovers were introduced. Still to be incorporated 
are dynamic allocation of gravity anomalies, surface densities. 
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and the new parameterization geoid heights. In the general area 
of modeling, the Jacchia 77 model has been implemented as well as 
a revised, solid-earth tidal model and general relativity model. 
The program now integrates in coordinate or ephemeris time 
instead of atomic time so that relativity can be properly 
modeled. The program has been rewritten for a vector or parallel 
processor . 

A gravity model has recently been developed using existing 
data to improve Seasat orbits for more accurate reduction of the 
spacecraft altimetry. Seasat orbits are now accurate to about 80 
cm in radial position. Further gravity model development is 
currently being done to improve Lageos orbits for more accurate 
baseline estimation. 

Long arc orbital solutions of Lageos laser ranging data for 
the May 1976 to January 1981 (Figure 13) interval revealed 
effects caused by ocean and solid earth tide model error and 
particular spherical harmonic coefficients of the gravity field. 

The estimation of GM, the product of the gravitational 
constant and the earth's mass, from monthly orbital arcs of 
Lageos, has continued through 1979. An error in the application 
of the 1976 data that was found earlier this year has been 
corrected and a complete set of 44 GM estimations has been 
obtained for the period May 1976 through December 1979. These 44 
values give GM = 398600.44 + .01 km^sec"^ compared with 39800.461 
+ 0.026 from the most recent lunar laser and Doppler results. 

The semimajor axis of Lageos' orbit is decreasing at the rate 
of 1.1 mm day" . The cause of this orbital decay appears to be a 
combination of charged particle and neutral particle drag. 
Calculations based on the charged and neutral particle density at 
the Lageos altitude plus laboratory charge drag experiments give 
the same order of magnitude as the observed orbital decay. Five 
other mechanisms (resonance, gravitational radiation, the 
Poynting-Robertson effect, transfer of spin to orbital angular 
momentum, and drag from near-earth dust) which can decrease the 
semimajor axis fail by at least an order of magnitude to explain 
the observed rate. Three other investigated mechanisms (the 
Yarkovsky effect, the Schact effect and terrestrial radiation 
pressure) also fail to explain a secular decrease in the 
semimajor axis. 

The passive geodetic satellite Starlette, which is in an 
orbit with an altitude of about 1000 km and an inclination of 
50°, is complementary to Lageos for precision geodynamic studies. 
Because of the lower orbital altitude, gravity modeling error 
causes ephemeris error of a few meters. Four years of precision 
laser tracking data (1975-1978) have been used to adjust the 
(36,36) GEM-lOB model plus additional terms out to degree 48. 
Tracking station coordinates were also adjusted in this solution. 
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A global 1° x 1° mean anomaly data tape used for geophysical 
studies and geoid computations, has been continuously updated and 
used in solutions combining satellite and gravimetric data. The 
latest update was made in October 1979, where 41973 anomalies 
were given. Of these values, 6912 anomalies were tentatively 
identified as being geophysically predicted. 

During the past year research was conducted to develop more 
accurate methods for the calculation of geoid undulations from 
potential coefficients and gravity anomaly data. Tests show that 
results are improved by a factor of two using a modification of 
the Molodensky truncation functions. 

The gravity field of the central Pacific area has been 
measured using satellite-to-satellite tracking between Geos-3 at 
an altitude of about 840 km and ATS-6 at six earth radii in a 
geosynchronous orbit. The measurements are about 70 km apart 
along-track. A low degree and order (less than 12) gravity field 
was assumed as the reference field in computing the orbits. 

Using about 50 tracks, these point values were contoured to give 
a gravity map at an altitude of 840 km. This field shows a 
dominant wavelength of about 2000 km, often with positive 
anomalies associated with residual depth anomalies. Significant 
new anomalies appear near the east Pacific rise trending in the 
ridge direction and continuing through North America. An 
additional group of about 50 new passes has been similarly 
reduced to produce a map which correlates closely with the newest 
GEM model, the Geos geoid, and with the Seasat geoid. 

The global sets of Geos-3 and Seasat altimeter data have been 
used to produce global 1° x 1° mean sea surface map (Figure 14) . 
Comparisons have been made between these surfaces and geoids 
computed from GEM gravity models. The rms difference between the 
Geos-3 and Seasat mean sea surfaces is 1.1 meters, indicating 
that sub-meter accuracy has been achieved for both surfaces. 
Adjustments of the origin of the reference ellipsoid for these 
surfaces indicated differences generally less than 50 cm. 


The Magsat spacecraft successfully acquired high-accuracy 
vector data from November 2, 1979, through June 11, 1980. Data 
analysis being performed includes modeling of the earth's core 
field and study of inferred core fluid properties, determination 
of characteristics of upper mantle conductivity, modeling of 
crustal features causing magnetic anomalies, and the study of 
fields from ionospheric and magnetospheric currents. 

An initial selection of Magsat data has been used to derive a 
series of magnetic field models. In addition to representing the 
core field well at the Magsat epoch, these models have been used 
to study the spectrum of the field to aid in separating core and 
crustal fields (Figure 15) . To model the field over longer time 
spans, and to provide a larger temporal data base for predicting 
the field, Magsat data has combined with pre-Magsat data from 
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both satellite and surface measurements to derive a 
time-dependent representation for the period 1960-1980. Surface 
data used in such models often has a very large contribution from 
local crustal sources which has not previously been accounted 
for. Because of the combination of satellites and surface data, 
the latest models include a solution representing the anomaly 
fields at observatories. Use of this technique improves not only 
the field representation but also the prediction ability of the 
resulting models. 

Preliminary examination of Magsat vector data indicates that 
it has the accuracy necessary to detect crustal magnetic anomaly 
fields. The vector data will be used to infer the presence or 
absence of a remanent component in the anomaly field. For 
example, a preliminary model of the Bangui anomaly indicates a 
good probability that it is entirely induced with no remanent 
component . 

The spherical harmonic spatial power spectrum for the 
geomagnetic field has been related to the sum of the squares of 
the coefficients of a spherical harmonic expansion of the scalar 
potential of the field, where the summation, for a particular 
degree, is over all orders m. The function is invariant under a 
rotation of the coordinate system, and is basic to a statistical 
description of the geomagnetic field. 

A statistical model for the crustal and core geomagnetic 
fields has been postulated, and equations have been derived 
giving the shape of the spherical harmonic spatial power spectrum 
for both the core and crustal geomagnetic fields. By fitting 
these equations to available data, equations have been obtained 
for the numerical values of the spatial power spectra for both 
the core and crustal geomagnetic field. The equations are not 
strongly model-dependent. The two spectra are approximately 
equal for n = 15. 

Equations have been derived relating the spherical harmonic 
spatial power spectrum to the average great circle Fourier power 
spectra of the magnetic field components. These equations have 
been used to estimate the crustal spherical harmonic spatial 
power structrum from data reported by Aldredge et al. in 1963 for 
an approximate great circle track flown by aircraft. It should 
be possible to use Magsat great circle data to improve this 
estimate without the need of doing a spherical harmonic analysis 
to high crustal frequencies. The equations can also be used to 
test for noise in the vector Magsat data. 


5. Advanced Geodynamics Instruments and Studies 

Gravsat results were simulated (one micron/sec tracking 
between two drag-free satellites at 200 km altitude, with 3° 
separation between them) show that harmonic coefficients (36 x 
36) can be obtained with an accuracy of 2 to 3 orders of 
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magnitude improvement over our present knowledge. The Gravsat 
harmonic solution was truncated at degree and order 36, showing 
increasing accuracy for the higher degree terms. However, 
results also show that to achieve this benefit the harmonic model 
must be carried out to very high degree (180 or higher) which can 
be accomplished only with the more advanced computers. The main 
benefit of the high accuracy in the low degree and order terms is 
that satellite position due to gravity modeling will be known to 
better than one centimeter. This is specially beneficial for 
estimating baselines from Lageos orbits and for estimating radial 
position on spacecraft orbits with altimetry application. 

An error analysis computer program was developed to compute 
the accuracy of mean anomaly and mean geoid undulations that 
would be expected from various types of Gravsat missions. The 
program assumes a polar orbit with global data coverage. The 
theory used is based on errors in the frequency domain and 
considerations of truncation effects. 

The possibility of using laser doppler tracking between the 
two Gravsat spacecraft has been studied; the approach would be to 
have a low-power helium-neon laser on each satellite, and to 
stabilize the lasers against temperature-controlled Fabry-Perot 
reference cavities. Telescopes 10 cm or less in aperture would 
be used. Preliminary results indicate that measurement 
accuracies of the order of 0.01 micron/sec may be obtained for 
periods up to at least 100 seconds, so that two orders of 
magnitude improvement could be made at shorter horizontal 
wavelengths, provided that comparable accuracy improvements could 
be made for measuring the along-track position of the optical 
systems with respect to the proof masses. 

Design studies have been conducted of modifying the 
Spaceborne Laser Ranging Concept to place the laser in an 
airplane rather than in a spacecraft (the passive retroref lectors 
are arrayed on the ground) . A system (Figure 16) is feasible 
which can range simultaneously to six ground-based 
retroref lectors in a 100 km by 200 km grid in a matter of hours. 
With range biases and single-shot rms errors of the order of one 
centimeter, such a system could determine baseline lengths of the 
order of 100 km to better than one centimeter. 

A second Lageos satellite would offer substantial benefits 
over the present situation. Studies show that station position 
errors can be reduced if the second satellite is placed in a 
posigrade orbit. A second satellite in the same orbit as Lageos 
(but 180° out of phase) does not offer great improvement in 
baseline length accuracies or station position, but does increase 
operational efficiency by reducing the time required to determine 
these quantities. 
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Small highly mobile systems based on use of the Global 
Positioning System satellites are being developed. The SERIES 
concept is to use standard VLBI techniques (Figure 17 ) , and it is 
expected that accuracies of the order of 1-2 cm can be achieved 
over baseline lengths of 2 to 200 km, with a few hours of on-site 
acquitition time. Quasar-referenced ARIES mobile facilities will 
be used to establish a sparse fundamental control grid, as a 
basis for making the GPS measurements traceable to the 
time-invariant quasar directions. Recent efforts have been 
directed toward detailed design and fabrication of a pair of 
proof -of -concept stations. A new ionospheric calibration system 
has been demonstrated. The SERIES stations contain lightweight 
1.5m diameter dish antennas mounted on trailers. It is expected 
that the two units will be tested in the field on baselines up to 
260 km in length, and their results compared with mobile laser 
ranging and quasar-based mobile VLBI stations. 

Geodynamic requirements for global reference frames have been 
studied, and problems have been identified with establishing 
services to maintain and publish reference frame information in 
the future, the reason being crustal movements which affect the 
station locations at the accuracies now being achieved. 

Monument stability is a fundamental limitation for many types 
of crustal movements measurements. An "optical anchor" has been 
installed at Pinon Flat Observatory in California; this device 
measures the movement of a laser strainmeter piers with respect 
to the material below. Significant noise level is observed with 
this instrument, which implies that similar equipment may be 
desirable in other observatories. 

A new method of determining polar motion is being studied; 
this uses a differencing procedure between quasi-simultaneous 
laser ranges. If successful, the method would reduce biases from 
systematic unmodeled orbital effects. 
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Appendix 1 

Geodynamics Investigations Funded Under Applications Notice ^ 1980 


Geopotential Fields 
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Paleo-Reconstructions . 

Langel, R. A., Goddard Space Flight Center: Magnetic Field 

Modeling, Including Secular Variation. 

Marsh, B. D., Johns Hopkins University: Interpretation and 

Delineation of the Earth's Gravity Field. 

Rapp, R. H., Ohio State University: Improvement of the 

Earth's Gravity Field from Terrestrial and Satellite Data. 


Global Earth Structure and Dynamics 


Daly, S. F., Jet Propulsion Laboratory: Variable Viscosity 

Convection. 

Dickman, S., SUNY Binghampton: Effects of the Ocean on Polar 

Motion. 

Farrell, W. E., Systems, Science & Software: Measurement of 

Earth Tides at NASA Tracking Stations. 

Felsentreger , T. L., Goddard Space Flight Center: Estimation 

of Solid Earth/Ocean Tidal Parameters from Satellite Tracking 
Data. 


Goodkind, J. M. , University of California at San Diego: 
Measurement of Vertical Crustal Motion with Superconducting 
Gravimeter . 

Haxby, W. F., Lamont-Doher ty Geological Observatory: The 

Effect of Subduction Zones on the Low-Order Harmonics of the 
Geoid . 

Kaula, W. M. , UCLA: Earth Dynamics Analyses. 
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Parmentier, E. M., Brown University; A Study of the Effects 
of Moving Continents; Continental Wakes. 

Schubert, G. , UCLA; Thermal Isostasy in the World's Oceans. 

Slade, M. A. , Jet Propulsion Laboratory; Post-Seismic Polar 
Motion Excitation - Numerical Studies. 

Smith, M. L., University of Colorado; Geophysical Studies of 
Nutation and Polar Motion. 

Williams, J. G. , Jet Propulsion Laboratory; Lunar Laser 
Ranging Modeling and Data Analysis. 

Yoder, C. F., Jet Propulsion Laboratory; Study of the Effect 
of Inner Core, Fluid Core, and Oceans on Earth Rotation, 

Nutation, and Polar Motion. 


Lithospheric Structure and Evolution 


Burke, K., SUNY Albany; Neotectonics of the Northern 
Caribbean Plate boundary Zone - A Study in Hispaniola. 

Dewey, J. F., SUNY Albany; Structural and Tectonic Evolution 
of the North Anatolian Fault; A Pilot Study in Synthetic 
Integration of Field Investigations and Landsat Interpretation in 
the Evolution of Intracontinental Plate Boundaries. 

Jones, G. , Texas A&M; Altimetry Data and the Deep Structure 
of Subduction Zones. 

Lowman, P., Goddard Space Flight Center; Late Cenozoic 
Tectonics of the Basin and Range Province. 

Molnar, P., MIT; Seismic Wave Velocity Structure in the 
Downgoing Slab & the Olivine-Spinel Phase Change. 

Roufosse, M. C., SAO; Global Study of the Time Evolution of 
the Lithosphere Using GEOS-3 and Seasat Altimeters. 

Watts, A. B. , Lamont-Doher ty Geological Observatory; The 
Relationship Between Gravity and Bathymetry in the Pacific Ocean. 
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Regional Crustal Deformation Modeling 


Berger, J., University of California at San Diego: 
Improvements in Benchmark Stability for Precision Tilt 
Measurements of Crustal Dynamics. 

Cohen, S., Goddard Space Flight Center: Models of Coseismic 

and Post-Seismic Deformation and Studies of Regional Seismicity. 

McAdoo, D. C., Goddard Space Flight Center: Dynamic 

Processes of Convergent Plate Margins. 

Melosh, H. J., SUNY Stony Brook: Finite Element 

Investigation of Lithosphere Tectonics. 

Molnar , P., MIT, Seismicity and Active Tectonics of the Andes 
and the Origin of the Altiplano. 

Mueller, I. I., Ohio State University: Advanced Studies for 

the Geodynamics Program. 

Nur, A., Stanford University: Models for Rupture Mechanics 

of Plate Boundaries and Crustal Deformation. 

Solomon, S. C., MIT, Application of Laser Ranging and VLBI 
Data to a Study of Plate Tectonic Driving Forces. 

Toksttz, M. N., MIT: Geophysical Study of the Structure and 

Processes of the Continental Convergence Zones Alpine - Himalayan 
Belt. 

Turcotte, D. L. , Cornell University: Finite Element Studies 

of the Surface Strain Field Adjacent to the San Andreas Fault. 

Whitcomb, J. H., CIRES (University of Colorado): Modeling of 

Tectonophysical Distortion from Measurements of Long-Baseline 
Geodetic Data and Other Geophysical Parameters. 


Advanced Studies 


Bender, P. L., University of Colorado: GPS Satellite Orbit 

Determination Using the Reconstructed Carrier Phase Method for 
Tracking . 

Counselman, C. C., MIT: Very Long Baseline Interferometric 

Geodesy with GPS Satellites. 

Resch, G. M. , Jet Propulsion Laboratory: Inherent Accuracy 

of the Water Vapor Radiometry Technique. 
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Glossary of Acronyms and Abbreviations 


AEM-C 

AFGL 

ALSEP 

AN 

AO 

ARIES 

ASTP 

ATS 

BIH 

GDP 

CIRES 

CSTG 

DMA 

DOC 

DOD 

DSN 

EDIS 

ESA 

FY 

GEOS 

GPS 

GSFC 

GUWG 

HO 

I AG 

lAU 

ICSU 

ICL 

ILP 

lUGG 

lUGS 

JPL 

Lageos 

LASSO 

LLR 

LURE 

Magsat 

Mark-III 

MERIT 

MIT 

MLRS 

Moblas 

MSFC 


Third Applications Explorer Mission (Magsat) 

Air Force Geophysical Laboratory 
Apollo Lunar Surface Experiment Package 
Applications Notice 
Announcement of Opportunity 

Astronomical Radio Interferometric Earth Surveying 

Apollo-Soyuz Test Project 

Applications Technology Satellite 

Bureau International de I'Heure 

Crustal Dynamics Project 

Cooperative Institute for Research in 

Environmental Sciences (University of Colorado) 
Commission for International Coordination of 

Space Techniques for Geodesy and Geodynamics 
Defense Mapping Agency 
Department of Commerce 
Department of Defense 
Deep Space Network 

Environmental Data Information Service 
European Space Agency 
Fiscal Year 

Geodynamic Experimental Ocean Satellite 
Global Positioning System 
Goddard Space Flight Center 
Gravsat User Working Group 
Haystack Observatory 
International Association of Geodesy 
International Astronomical Union 
International Council of Scientific Unions 
Inter-Union Commission on the Lithospehre 
International Lithosphere Program 
International Union of Geodesy and Geophysics 
International Union of Geological Sciences 
Jet Propulsion Laboratory 
Laser Geodynamics Satellite 

Laser Synchronization (of atomic clocks) from 
Synchronous Orbit 
Lunar laser ranging 
Lunar Ranging Experiment 
Magnetic Field Mapping Satellite 
Advanced VLBI data system 
Monitoring Earth Rotation and Inter- 
comparison of Techniques 
Massachusetts Institute of Technology 
McDonald Laser Ranging System 
Mobile Laser System 
Marshall Space Flight Center 
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NAS/NRC 

NASA 

Natmap 

N3S 

NGS 

NEROC 

NGSP 

NOAA 

NRAO 

NSF 

NSSDC 

NSWC 

OSTA 

OVRO 

Polaris 

PPME 

Ramlas 

RTD 

RTOP 

SAFE 

SAG 

SERIES 

Stalas 

STDN 

TLRS 

TPM 

USGS 

USNO 

VLB I 


National Academy of Sciences - National 
Research Council 

National Aeronautics and Space Administration 
Division of National Mapping (Australia) 

National Bureau of Standards 
National Geodetic Survey 

Northeastern Radio Observatory Corporation 
National Geodetic Satellite Program 
National Oceanic and Atmospheric 
Administration 

National Radio Astronomy Observatory 

National Science Foundation 

National Space Science Data Center 

Naval Surface Weapons Center 

Office of Space and Terrestrial Applications 

Owens Valley Radio Observatory 

Polar Motion Analysis by Radio Interferometric Systems 
Pacific Plate Motion Experiment 
Range Measurement by Laser System 
Research and Technology Development 
Research and Technology Operating Plan 
San Andreas Fault Experiment 
Smithsonian Astrophysical Observatory 
Satellite Emission Radio Interferometric 
Earth Surveying 
Stationary Laser System 
Space Tracking and Data Network 
Transportable Laser Ranging System 
Tectonic Plate Motion 
United States Geological Survey 
United States Naval Observatory 
Very Long Baseline Interferometry 
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Appendix 3 


Publications of the Geodynamics Program, 1976-1981 


Adams, J., Cornell University: Active Tilting of the 

Midcontinent: Geodetic and Geomorphic Evidence. EOS Trans. 

Amer. Geophys. Union 60, 310, 1979. 

Adams, J., and Reilinger, R.E., Cornell University: Time 
Behavior of Vertical Crustal Movements Measured by Releveling in 
North America. Proceedings of the International Symposium on 
Redefinition of the North American V ertical Geodetic Networks, G. 
Lachapelle, ed., Ottawa, 327-342, 1980. 

Adams, J., and Reilinger, R.E., Cornell University: Active 

Tilting of the Oregon and Washington Coastal Ranges. EOS Trans. 
Amer. Geophys. Union 61, 371, 1980. 

Afonso, G., Centre d'Etudes et de Recherches, Geodynamiques 
et As tronomiques : Role majeur de la trainee electrique dans la 

decroissance seculaire du demi-grand axe du satellite Lageos. 
These, Univ. Paris, 1980. 

Afonso, G. , Centre d'Etudes et de Recherches, Geodynamiques 
et As tronomiques: Effect of Atmospheric and Charged-Particle 

Drag on the Lageos Satellite. C.R. Acad. Sci. Paris 290, 1980. 

Alexander, S.S., Lee, B., and Petrus, C., The Pennsylvania 
State University: Relation of Earthquake Focal Mechanisms to 

Tectonic Structure Features Observed on Satellite Imagery. EOS, 
Trans Amer. Geophys. Union 61, 290, 1980. 

Banerjee, B., The Johns Hopkins University: Modeling the 

Gravity Field of the Earth by Point Masses. EOS Trans. Amer. 
Geophys. Union 61, 371, 1980. ' 

Biadgelgne, A., The Pennsylvania State University: Third 

Derivative Analysis of Gravity Anomalies From Prismatic Sources. 
MS Thesis, The Pennsylvania State University, November 1979. 

Bock, Y., Mueller, 1. 1., and Pavlis, E., Ohio State 
University: On the VLBI-Satellite Laser Ranging 'Iron Triangle' 

Intercomparison Experiment. Radio Interferometry Techniques for 
Geodesy, NASA Conference Publ. 2115, 435-445, 1980. 

Bock, Y., Ohio State University: A VLBI Variance— Covariance 

Analysis Interactive Computer Program. Report 298, Dept. 

Geodetic Science, Ohio State University, 1980. 
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Bowin, Carl 0., Woods Hole Oceanographic Institution: Mass 

Anomaly Estimation From Combined Inversion of Gravity and Geoid 
Data. EOS Trans. Amer. Geophys. Union, 1981. 

Bowin, Carl 0. , Woods Hole Oceanographic Institution; Upper 
Mantle Heterogeneity; Comparison of Region South of Australia 
With Philippine Basin. EOS Trans. Amer. Geophys. Union, 1981. 

Bowman, P.L., Purdue University; Correlation of Gravity 
and Magnetic Data Over North America. MS Thesis, Purdue 
University, 1978. 

Bowman, P.L., Braile, L.W., Chandler, V.W. , Hinze, W.J., 

Luca, A.J. , and von Frese; Magnetic and Gravity Anomaly 
Correlation and its Application to Satellite Data. NASA 
TM-79702, 1979. 

Brammer, R.F., Sailor, R.V. , and Lazarewicz, A.R. , The 
Analytical Sciences Corporation (TASC) ; Magsat Investigation of 
Lithospheric Magnetic Anomalies in the Eastern Indian Ocean. EOS 
Trans. Amer. Geophys. Union, May 1981. 

Brov/n, L.D., Miesen, D.L., Reilinger, R.E. , and Jurkowski, 
G.A., Cornell University; Geodetic Leveling and Crustal Movement 
in the U.S., Part I, Topography and Vertical Motion. EOS Trans. 
Amer. Geophys. Union, changes and corrections for Spring Annual 
Meeting, page 4, 1980. 

Carter, W.E., Rogers, A.E., Counselman, C.C., and Shapiro, 

I. I., National Geodetic Survey, National Ocean Survey, NOAA: 
Comparison of Geodetic and Radio Interferometric Measurements 
of the Haystack-Westford Base Line Vector. J. Geophys. Res. 85, 
1684-2687, 1980. 

Chaffin, David L., The Pennsylvania State University; 
Implications of Regional Gravity and Magnetic Data for 
Structures Beneath Western Pennsylvania. MS Thesis, The 
Pennsylvania State University, June 1981. 

Chandler, V.W., Bowman, P.L., Hinze, W.J., and O'Hara, N.W. , 
Purdue University; Long-Wavelength Gravity and Magnetic 
Anomalies in the Lake Superior Region. Geol. Soc. Amer. Special 
Paper on Geology and Tectonics of the Lake Superior Basin, in 
press, 1981. 

Chandler, V.W., Hinze, W.J., and Braile, L.W., Purdue 
Uniyersity: A Preliminary Study of Combined Gravity and Magnetic 
Analysis. EOS Trans. Amer. Geophys. Union 56, 597, 1975. 

Chandler, V.W. , Hinze, W.J., and O'Hara, N.W., Purdue 
Uniyersity; Regional Analysis of Magnetic and Gravity Anomaly 
Data in the Great Lakes Region, North America. EOS Trans. Amer. 
Geophys. Union 59, 378, 1978. 
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Chandler, V.W., Hinze, W.J., and Braile, L.W., Purdue 
University: Correlation of Gravity and Magnetic Data Over the 

Central United States. EOS Trans. Amer. Geophys. Union 59, 226, 
1978. 

Chandler, V.W. , Koski, J.S., Hinze, W.J., and Braile, L.W., 
Purdue University: Analysis of Multi-source Gravity and Magnetic 

Anomaly Data Sets by Moving-Window Applications of Poisson's 
Theorem. Geophysics 46, 30-39, 1981. 

Chandler, V.W. , Purdue University: Correlation of Gravity 

and Magnetic Data Over the Great Lakes Region, North America. 
Ph.D. Thesis, Purdue University, 1977. 

Chen, C.H., Chen, W.P., and Molnar, P., Massachusetts 
Institute of Technology: The Uppermost Mantle P Wave Velocities 
Beneath Turkey and Iran. Geophys. Res. Lett 7, 77-80, 1980. 

Chen, W.P., and Molnar, P., Massachusetts Institute of 
Technology: Constraints on the Seismic Wave Velocity Structure 

Beneath the Tibetan Plateau and Their Tectonic Implications. J . 
Geophys. Res., in press, 1981. 

Chen, W.P., and Molnar, P., Massachusetts Institute of 
Technology: Constraints on the Seismic Wave Velocity Structure 

Beneath the Tibetan Plateau. Proceedings of the Symposium on the 
Tibetan Plateau, Beijing, in press, 1981. 

Cheney, R.E., and Marsh, J.G., Goddard Space Flight Center: 
Seasat Altimetry Observations of Dynamic Ocean Currents in the 
Gulf Stream Region. J. Geophys. Res., submitted, 1980. 

Cheney, R.E., and Marsh, J.G., Goddard Space Flight Center: 

A Preliminary Look at Seasat Altimetry in the Gulf Stream Region. 
Gulf stream 5, 1979. 

Chi, S.C., Reilinger, R.E., Brown, L.D., and Oliver, J.E., 
Cornell University: Leveling Circuit Closures and Vertical 

Crustal Movements. J. Geophys. Res. 85, 1469-1474, 1980. 

Chi, S.C., Reilinger, R.E., Brown, L.D., and Jurkowski, G.A., 
Cornell University: Geodetic Leveling and Crustal Movement in the 
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Geophys. Union 60, 315, 1980. 

Chr istodoulidis, D.C., and Smith, D.E., Goddard Space Flight 
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Amer. Geophys. Union 61, 1980. 
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Figure Captions 


1. Mobile VLBI station MV-3. 

2. TLRS-1 during checkout at the University of Texas at Austin. 

3. TLRS-2 in the laboratory at the Goddard Space Flight Center. 

4. VLBI/laser intercomparison sites. Figures on each baseline 
are deviations (laser minus VLBI) between baseline measurments, 
in centimeters. 

5. Sites planned for occupation to measure regional deformation 
in North America, 1980-86. 

6. Sites planned for occupation to measure crustal deformation 
in California. 

7. Artist conception of Magsat spacecraft in orbit. Magneto- 
meters were mounted on the extension boom. 

8. Gravsat satellite concept. Satellite-to-satellite tracking 
antennas protrude above and below the spacecraft. The second 
spacecraft is in the background. TDRSS antennas are the strips 
across the nose and tail of the spacecraft. 

9. Global atmospheric angular momentum compared to length- 
of-day observations, January 1976 through April 1980. 

10. Polar motion: comparison between VLBI and Lageos laser 

ranging measurements during 1980 MERIT observing campaign. (a) 
x-component of pole position; (b) y-component of pole position. 

11. Baseline length changes, Quincy to Otay Mountain 
(California) , from San Andreas Fault Experiment, 1972-1979. 

12. Baseline length residuals. Haystack Observatory 
(Massachusetts) to Owens Valley Radio Observatory (California) 
September 1976 through September 1980. 

13. Polar motion and earth rotation based on Lageos 
observations. May 1976 through December 1979. 

14. Mean sea surface topography based on Seasat altimeter data. 

15. Magsat magnetic anomaly map. 

16. Airborne Laser Ranging System block diagram. 

17. Prototype SERIES unit at the Jet Propulsion Laboratory. 
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